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^3 Abstract. One of the most exciting potential sources of gravitational waves for the 

Laser Interferometer Space Antenna (LISA) are the inspirals of approximately solar 

00 mass compact objects into massive black holes in the centres of galaxies — extreme 

mass ratio inspirals (EMRIs). LISA should observe between a few tens and a few 

(J hundred EMRIs over the mission lifetime, mostly at low redshifts (z < 1). Each 

observation will provide a measurement of the parameters of the host system to 
%—i unprecendented precision. LISA EMRI observations will thus offer a new and unique 

1 way to probe black holes at low redshift. In this article we provide a description of 

the population of EMRI events that LISA is likely to observe, and describe how the 
numbers of events vary with changes in the underlying assumptions about the black 

qq hole population. We also provide fitting functions that characterise LISA's ability to 

00 detect EMRIs and which will allow LISA event rates to be computed for arbitrary 

population models. We finish with a discussion of an ongoing programme that will 
use these results to assess what constraints LISA observations could place on galaxy 
evolution models. 
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1. Introduction 

The inspiral of a stellar mass compact object — a black hole (BH), neutron star (NS) or 
white dwarf (WD) — into a massive black hole (MBH) with mass in the range 1O 4 M - 
1O 7 M in the centre of a galaxy will generate gravitational waves at frequencies to which 
the planned space based gravitational wave detector, LISA [1] will be sensitive. These 
"extreme mass ratio inspiral" (EMRI) sources are of particular interest to relativists 
as the emitted gravitational waves encode detailed information about the spacetime 
structure close to the massive central object and hence can be used to test whether these 
objects are indeed the Kerr black holes predicted by relativity (see [2] and references 
therein). However, these sources are also of great interest for astrophysics. A LISA 
observation of an EMRI event can determine the parameters of the system to very high 
precision — the mass of the inspiraling object and the mass and spin of the central black 
hole are typically recovered to accuracies of a fraction of a per cent |3j. Preliminary 
estimates have indicated that LISA may see as many as a few thousand EMRI events over 
the mission lifetime [4j . This set of observed EMRI events will provide us with detailed 
information about the properties of black holes in the relatively nearby Universe. 

To date, EMRI event rate calculations have been done only crudely and for specific 
models of the black hole population. In addition, there is no simple prescription in 
the literature to go from an astrophysical population model to a realisation of the LISA 
event distribution. In this paper, we will describe some of the properties of the likely set 
of LISA EMRI events more carefully. In order to do so, we have computed the distance 
to which LISA will see different types of events and have folded these ranges into a 
simple population model to calculate the expected number and distribution of observed 
events. In the course of this work, we have derived fits which give the detectable depth 
of an EMRI event as a function of its parameters and we also provide these here. We 
present EMRI event distributions for a simple set of MBH population models in which 
there is no evolution over the redshift range to which LISA is sensitive (0 < z < 1.5). 
These models illustrate the potential for astrophysics that LISA EMRI observations 
offer, while the fits will enable other researchers to examine alternative galactic black 
hole population models and readily convert these into LISA event distributions. 

The paper is organised as follows. In Section [2] we describe the waveform model 
we have used to compute LISA signal-to-noise ratios and provide fits to the "observable 
lifetime" of events, which characterises the likelihood that LISA could see them. In 
Section [3] we use these results to compute event rates assuming a simple non-evolving 
MBH population. We compute the distribution of LISA events, and describe how the 
number of events and the minimum and maximum likely redshift of observed events vary 
with the inspiral rate and the spin of the MBH. Finally, in Section [4] we discuss how 
we are currently using these results to investigate LISA's ability to distinguish between 
different galaxy evolution models, as well as describing some of the limitations of the 
current work and how we plan to address these in the future. 
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2. LISA sensitivity to EMRIs 

To compute the population of LISA EMRI events, we need a prescription to determine 
when an event is detectable. Typically, we expect an event to be detected if the 
matched filtering signal-to-noise ratio (SNR) exceeds some predetermined threshold, 
Pthresh- The first estimate of the SNR threshold that will be needed was pthresh ~ 35, 
which was computed assuming that the data analysis would be carried out using a 
hierarchical semi-coherent matched filtering algorithm [4J. This is somewhat higher than 
the threshold estimated to be needed for a fully coherent search, pthresh ~ 15, if such 
a search was computationally feasible. Several alternative approaches to EMRI data 
analysis have been proposed subsequently, including time-frequency analyses [5j EJ [U [8] 
and Metropolis-Hastings Monte Carlo (MHMC) searches [HJ [TUl ITT] . Time-frequency 
analyses will require higher SNRs for detections. MHMC searches may be able to get 
closer to the fully-coherent limit of pthresh ~ 15, but this has not yet been demonstrated 
for isolated sources let alone multiple overlapping sources. In this work we will take 
Pthresh = 30 and assume that the LISA mission lasts 5 years, with no failures on the 
satellite so that the two-independent low-frequency channels [12] are available for the 
entire mission. This threshold may be an overestimate of what will be required, but the 
other assumptions are optimistic, as it is unlikely that LISA will have full functionality 
for 5 years. Results will be based on this LISA configuration ('5yr, 2d') unless explicitly 
stated otherwise, but this should be regarded as optimistic. Given the uncertainty in 
data analysis practicalities and on the lifetime of the mission, we will also present results 
from a pessimistic case in which we assume the mission lasts only two years and there 
is a partial failure on the satellite which means only one of the low-frequency channels 
is available ( £ 2yr, Id'). 

If EMRI events were short- duration, the event rate would be determined by 
computing the distance at which the SNR equals pthresh and then multiplying the rate 
per unit volume by the volume contained by that distance. However, EMRIs are long- 
lived, and SNR can be accumulated for as much of the inspiral as coincides with the 
LISA observation. For a given system at a given distance it is possible to compute the 
SNR as a function of the time remaining until plunge, T pi . If the LISA mission lasts T L isa 
years, then for t p \ < Tlisa, this SNR is accumulated over r p i years, while for r p i > Tlisa, 
this SNR is accumulated over the Tlisa years from r p \ = r to t p i = t — Tlisa- As 
a general rule, the SNR will initially increase as r p i decreases and then decrease as t p i 
approaches zero. So, there will be a largest (earliest) r p \ for which the SNR exceeds 
Pthresh, T eari y say, and a smallest (latest) r ph r late say. If the LISA satellite turns on when 
a particular system is at any point in the range n a te < T pi < Nearly, that system will be 
seen by LISA. We can thus define an observable lifetime for EMRI systems, r, by 

r(A) = r oar i y (A) - r latc (A) (1) 

where A denotes the system parameters. If EMRIs plunge at a rate TZ per year in a 
particular galaxy, then tTZ gives the expected number of events that LISA will observe 
from that galaxy, after appropriate averaging of 1Z over the other system parameters, A. 
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2.1. SNR Calculation 

To compute SNRs, we use the fluxes presented in Finn and Thorne |13| . They computed 
the gravitational wave emission from circular and equatorial EMRIs via solution of the 
Teukolsky equation, i.e., accurate to arbitrary orders in velocity, but only to leading 
order in mass-ratio. These are the most accurate fluxes presently available for a range of 
central black hole spins. The restriction to circular and equatorial orbits is unphysical 
in that most mechanisms predict EMRI orbits will be of moderate eccentricity and 
inclined to the black hole equatorial plane when they enter the LISA band [2] . However, 
consideration of circular-equatorial orbits is a necessary first step and should indicate 
what the dependence of the LISA sensitivity on the central black hole mass and spin 
will be like, the investigation of which is the primary goal of this paper. Eccentricity 
and inclination introduce complexity into the waveform which will impact these results. 
The extension of this work to eccentric and inclined orbits is a necessary future project, 
which could be undertaken using approximate EMRI waveform models [31 H3] . 

The SNR, p, of an EMRI consisting of a compact object of mass p falling into a 
massive black hole of mass M and dimensionless spin a = S/M 2 at a redshift z is given 

by m 



P 2 
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In this, S% A denotes the sky and orientation averaged spectral density of the detector, 
di(z) is the luminosity distance to redshift z and = M(l + z)Vt = l/(a + (r/M) 3 / 2 ) is 
the dimensionless angular frequency of the orbit when the object is at Boyer-Lindquist 
radius r. The 7i Cjm 's are relativistic correction factors to the post-Newtonian waveform 
amplitudes that may be computed from tables in [13]. The summation over waveform 
harmonics stops at m = 4, since no higher harmonics are tabulated in [13]. However, 
the SNR contributed by the m = 4 harmonic is already a small fraction of the total, 
and so the inclusion of harmonics with m > 5 will make only a small correction to 
these results. Finn and Thorne tabulate results only as far as r/r isco = 10. Beyond 
that radius, we use extrapolations based on polynomials in \j\fr with the constraint 
that Ti^m — >• 1 as r — > oo. We have checked that our results are not sensitive to the 
exact way in which the extrapolation is done. In this analysis, we used S^ A , including 
confusion from white-dwarf binaries, taken from [3]. For the '5yr, 2d' case, we multiplied 
the SNR by a factor of \/2 to account for the fact that LISA can be thought to consist 
of two independent detectors at low-frequency. We computed d^z) using a standard 
A-CDM cosmology with H = 71km s" 1 Mpc" 1 , tt M = 0.27 and fi A = 0.73. These 
SNRs use a low-frequency approximation to the LISA response, and so we will slightly 
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underestimate the SNRs for systems with smaller central black hole masses, M, since 
at higher frequencies the LISA response gains a third independent channel. 

Expression ^ allows us to compute the observable lifetime for EMRIs as a function 
of /i, M, a and z. We note that we are computing the observable lifetime of a sky-and- 
orientation-averaged source, rather than the sky-and-orientation-averaged observable 
lifetime of a source, which is what we actually need, but we hope these will be 
comparable. 



2.2. Observable lifetimes 

In Figure [I] we show contours of constant observable lifetime on the (M, z) plane for 
black hole inspirals with [i = 1OM . Note that here and elsewhere we quote masses and 
observable lifetimes as measured at the source as these are of most use when computing 
event rates, since the rate of EMRIs per black hole, 71, is usually quoted in the source 
frame. The apparent mass and observable lifetime will be M(l + z) and r(l + z). We 
see that LISA is most sensitive to central black holes with mass around 1O 6 M , and 
can see many systems out to z ~ 1, with prograde inspirals into rapidly spinning MBHs 
being visible out to z ~ 2. For rapidly spinning black holes, the distance sensitivity is 
increased and shifted to higher intrinsic MBH masses. This follows from the fact that 
the emission comes out at higher frequencies for a given black hole mass, and hence the 
frequency at the floor of the LISA sensitivity curve corresponds to larger M. 

The lifetime r, for a fixed redshift z, can be well approximated by a trapezium 
function of log(M), with break-points at log(M) = x m i n , x\, x 2 , s max and lifetimes of 
r = 0,|/i,2/2) at those points; i.e., 

log(M) < x min 

y 1 (log(M) - x min )/ (xt - x mhl ) x min < log(M) < x x 

t < yi + (2/2 - Vi) (log(M) - x 1 )/(x 2 - zi) x 1 < log(M) < x 2 .(3) 

y 2 (1 - (log(M) - x 2 )/(x max - x 2 )) x 2 < log(M) < x max 

x max < log(M) 

This is illustrated in Figure [2] for a few cases. For low z and high spin in particular, the 
observable lifetime shows a lot more features which arise due to the shape of the LISA 
noise curve. However, we ignore these in the interest of having a simple fitting function. 
We will see in Figure [3] that the fits reproduce the event distribution well. We have also 
obtained fits to the six functions x m i n , xi, yi,x 2 , y 2 , x max as a function of z for the different 
central black hole spins quoted in [13] . The fits are given in Table [l] for a = 0,0.9,-0.9. 
Corresponding fits for the pessimistic detector configuration are given in Table [2j Note 
that here and elsewhere, an inspiral with a < denotes a retrograde inspiral into a 
black hole of spin \a\. For other values of the spin, a reasonable approximation may 
be obtained by linear extrapolation of the lifetimes given here. The functional forms of 
these fits were chosen arbitrarily, and these functions could probably be fit equally well 
using a different prescription. Given the other uncertainties about the EMRI population, 
these fits are probably more accurate than required in practice. However, we have opted 
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Figure 1. Contours of constant observable lifetime for BH inspirals. The upper plot 
shows contours with t = lyr for a range of central black hole spins, while the lower 
plot shows contours for prograde orbits into central black holes with spin a — 0.9 for 
a range of values of r. 

to provide an accurate description of the instrumental response (which is less uncertain 
than the astrophysics) and encourage readers to simplify these fits or not according to 
the application. We hope that these fits will be useful for computations of LISA EMRI 
event distributions for particular MBH population models. 

The fits in Tables [Ijj2] apply to inspirals of black holes (// = 1OM ). We provide only 
these results as such events will dominate the LISA event rate. The observable lifetime 
is affected by the shape of the LISA noise curve and the imposed SNR threshold, so 
the scaling with /x is non-trivial. However, for lower values of \i the observable lifetime 
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Figure 2. Trapezium fits to the observable lifetime as a function of central black 
hole mass, for a variety of spins and source redshifts, as indicated in the legend, and 
with fi — 10 Mq. Lifetimes computed from the numerical observable lifetime data are 
indicated as "Raw" and those computed from the fit by "Fit". The values of a and 
z for this plot were chosen at random, but these curves are representative of all other 
choices. We see that the simple trapezium fit works well in all cases. 



tends to be longer at lower redshift, since the inspiral proceeds more slowly, but the 
observable lifetime falls off more quickly as the redshift increases, since the instantaneous 
amplitude is lower and the threshold SNR is reached at a much lower redshift. As 
an illustration, the observable lifetime (in years) for a = 0.9 and M = 1O 6 M 
at z = {0.01,0.05,0.1,0.2,0.3,0.4,0.5} is {31,23,14,8.7,6.5,5.1,4.3} for y, = 1OM , 
{91,43,21,8.3,4,1,0} for y = 1.4M and {187,40,14,0.6,0,0,0} for fj, = O.6M . 
For M = 3 x 10 5 these numbers become {33,26,17,10,6.8,5,4} for y = 1OM , 
{111,42,16,4,0,0,0} for /i = 1.4M and {197,28,4.5,0,0,0,0} for // = O.6M . 



3. LISA EMRI event rates 



To compute the number and parameter distribution of events that LISA will see, we 
need a prescription for the comoving number density of MBHs, J\f(M, a, z), and the rate 
of EMR inspirals in any particular system, 1Z(M, a, z). The quantity J\f(M, a, z)dMda 
denotes the number of massive black holes per comoving volume with mass in the range 
M — > M + dM and spin in the range a — > a + da. As little is known about black 
hole spins, we assume that the MBH mass and spin distributions are independent and 
factorize J\f(M, a, z)dMda = (dn/dln M)(M, z) din M p(a, z)da, with normalisation 
J p(a,z)da = 1. For simplicity we further assume that there is no evolution in the 
properties of black holes with z, so that (dn/dln M)(M, z) = (dn/d In M)(M) and 
p(a, z) = p(a) only. This is probably reasonable for the range of redshifts of LISA 
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Table 1. Fitting functions for coefficients in the trapezium fit to the observable 
lifetime described by Eq. |3|. The fits are valid for z < 1, z < 1.45 and z < 0.7 for 
a = 0, 0.9, —0.9 respectively. For z outside this range, r = 0. 



EMRI events. The number of LISA events is then given by 

N L1SA = / / / n(M)r(M,a,z)-— -(M)p(a)-^ da din Mdz. (4) 

J z=0 Jm=m 1ov Ja=~\ d In M dz 

Here (dV c /dz) dz is the comoving volume in the redshift range z to z + dz, which we 
compute for the same A-CDM cosmology used previously. 

The mass function, dn/dln(M), can be derived from observed galaxy luminosity 
functions using the L — a and M — a relations. Using results from Aller & Richstone [T5] , 
one finds the function dn/dlnM is approximately flat and equal to 3.6 x lO^/i^Mpc -3 
for MBHs in the LISA mass range, M < 5 x 1O 6 M (NB h n = F /(71km s" 1 Mpc" 1 )). 
If we remove Sc-Sd galaxies (several of these are known to have black holes of much 
lower mass than that derived from the luminosity [IS]), this becomes approximately 
1.8 x lO^/i^Mpc -3 . This was the expression used to derive LISA EMRI event rates 
in [3]. However, there is great uncertainty in the black hole mass function in the LISA 
range, since there are very few measurements of the masses of quiescent black holes in 
that range [TTJUS]- We therefore adopt the prescription 
dn _ ( M 
dlnM ~ n ° \3x 10 6 M ( 

and will normalise to the case no = 0.002Mpc~ 3 and a = 0. 
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Table 2. As Table [T] but now for the pessimistic detector case, '2yr, Id'. The fits are 
valid for z < 0.55, z < 1 and z < 0.42 for a = 0,0.9, —0.9 respectively. For z outside 
this range, r = 0. 



We can adopt a similar ansatz for the intrinsic EMRI rate per galaxy, i.e., the 
frequency with which EMR inspirals start or end within a particular system 

/ M 

Here 7 denotes the type of EMRI — BH, NS or WD (which we take to have masses 
fi = 1OM , 1.4M , 0.6M Q respectively) - - and the subscript 'MW indicates we are 
normalising to a mass comparable to the Milky Way black hole. Hopman [T9] 
quotes rates and scalings of 7?.^ w = TZq = 400Gyr~ 1 , 7Gyr~ 1 , 20Gyr _1 and /3 7 = 
—0.15, —0.25, —0.25 for inspiraling black holes, neutron stars and white dwarfs 
respectively. The subscript '0' is just used to distinguish these fixed reference rates 
from our chosen 7£^ w . There is presently no observational prescription for the MBH 
spin distribution, p(a), so we use the simple ansatz p = (5 (a — a ) + 5 (a + a ))/2 for 
three cases a = 0, 0.5, 0.9. This assumes prograde and retrograde inspirals are equally 
likely in a system with given spin. 

Putting this together we can quote the number of events per species as 

N Njf 7 ^) f y ( a + P y ) ( ^ V (7) 

7 V ^0 / V0.002Mpc- 3 y v ; 
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The coefficients Nq and functions / 7 are tabulated in Table [3] for both the optimistic 
and pessimistic detector configurations. We normalise the rates per galactic black hole 
to the reference rates of Hopman [19] and normalise the / 7 's such that / 7 = 1 when 
a + /3 7 = -0.15, -0.25, -0.25 for BHs, NSs and WDs respectively. Note that in this 
simple prescription we cannot disentangle the effects of a and /3 7 on the event rate. We 
note also that we only include MBHs in the range 1O 4 M < M < 3 x 10 7 M Q in our event 
rate computation. The main conclusion from this table is that the event rate is very high 
and is totally dominated by BH inspirals, since these are visible much further than NS 
or WD inspirals, and are intrinsically more common in the reference model. Only if the 
WD or NS rates are enhanced by other mechanisms (e.g., tidal stripping of binaries [2]) 
or rise much more steeply toward low M will we see significant numbers of those events. 
Using the pessimistic detector assumptions, the rate decreases by approximately a factor 
of 10, although the rate scaling, / 7 , does not change significantly. The event rates 
in the optimistic detector case are quite large compared to previous EMRI event rate 
estimates jl]. This is primarily because earlier work took the rate scaling to be /3 7 = 3/8 
for all species, which suppresses the number of events with M < 3 x 1O 6 M and these 
are the dominant contributors to our event rate. With /3 7 = 3/8 and a = 0.9 we predict 
0.7, 1.2 and 650 events for WDs, NSs and BHs respectively. The number of WD events 
(~ 100) predicted in [I] was two orders of magnitude greater than the result we quote 
here. This arose because this previous work used a much higher WD rate of 5Myr _1 , 
which was based on older simulations. All of these numbers should be treated with 
caution, as there are very large uncertainties in the intrinsic rate of inspirals per galaxy, 
and the rates we have used could be off by as much as two orders of magnitude [2]. At a 
rate of 4 BH inspirals per Gyr per galaxy, we will see only 10-20 events over the mission 
lifetime, even with optimistic assumptions about the detector. 

In Figure [3] we show the distribution of events as a function of M and z for the 
BH events in our reference system with ao = 0.9. The distribution is peaked at a mass 
slightly below 1O 6 M , and events are mostly confined to the range 10 5 M Q < M < 
2 x 1O 6 M , although there is a significant tail toward smaller central black holes. The 
redshift distribution is peaked at z ~ 0.4 with a long tail out to redshifts z > 1. With 
pessimistic detector assumptions, the distribution is shifted to lower intrinsic masses 
and lower redshifts, and with a factor of ten fewer events overall. A useful quantity to 
compute is the value of z, z min , such that the expected number of events with z < z min is 
equal to 1. This is an estimate of the lowest redshift source that we might observe. The 
highest likely redshift of a source can be defined in a similar way. We tabulate these 
quantities for several different BH EMRI rates in Table |4| It is clear that we expect 
most of our events to be in the range 0.1 < z < 1, but if MBHs tend to have high spins 
and the intrinsic EMRI rate per MBH is also high, we could see an event with redshift 
as high as z = 1.5 or as low as z = 0.02. We note that the minimal likely redshift is 
independent of the central black hole spin. This is because nearby sources are visible for 
longer, and inspirals spend proportionally longer far from the central black hole. The 
nearest source detected is therefore likely to be observed early, i.e., long before plunge, 
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6 

[0.06] 


PYnf- fl fi9 — 9 1P,t 4- fl 47t 2 1 

[exp(-0.57 - 2.18s + 0.36a; 2 )] 


0.5 


BH 


[170] 


e.x.p^ — u.oy — z.o^x ^ u.u^a j 
[exp(-0.37 - 2.39a; + 0.60a; 2 )] 


NS 


1 4 
[0.1] 


p-xrif— 56 — 9 08r -1- 57r 2 — 081 t 3 ) 
[exp(-0.51 - 1.93a; + 0.46a; 2 - 0.044a; 3 )] 


WD 


n 75 
[0.07] 


pxnf — 5S — 1 Q8r -I- 54r 2 — 081 r 3 ) 
[exp(-0.49 - 1.86a; + 0.43a; 2 - 0.037a; 3 )] 


0.9 


BH 


1580 
[260] 


exDf— 30 — 1 85x + 82a; 2 — Q73x 3 ) 
[exp(-0.26 - 1.60a; + 0.77a; 2 - O.lOx 3 )] 


NS 


2.6 
[0.3] 


exp(-0.36 - 1.27a; + 0.70a; 2 - 0.14a; 3 ) 
[exp(-0.32 - 1.11s + 0.59a: 2 - 0.10a; 3 )] 


WD 


1.4 
[0.1] 


exp(-0.33 - 1.13a; + 0.66a; 2 - 0.14a; 3 ) 
[exp(-0.29 - 1.02a; + 0.54a; 2 - 0.091a; 3 )] 



Table 3. Normalisation factors and functions appearing in Eq. (J7J. The dependence 
on the rate scaling, a + /3 7 , was obtained by empirical fitting with simple functions 
in the range —1.5 < x < 1.5 and may not be trustworthy outside this range. The 
second line in each row, denoted by square brackets, is the corresponding result for the 
pessimistic '2yr, Id' configuration. 



Spin 





0.5 


0.9 


^(Gyr- 1 ) 


4 


40 


400 


4 


40 


400 


4 


40 


400 


5yr, 2d 
2yr, Id 


0.14 
0.30* 


0.05 
0.07 


0.02 
0.03 


0.14 
0.29* 


0.05 
0.07 


0.02 
0.03 


0.14 
0.27* 


0.05 
0.07 


0.02 
0.03 


5yr, 2d 
W 2yr, Id 


0.64 
0.16* 


0.76 
0.41 


0.84 
0.49 


0.74 
0.19* 


0.94 
0.49 


1.0 
0.59 


1.1 
0.41* 


1.3 
0.76 


1.4 
0.86 



Table 4. Most likely minimum and maximum redshift LISA events that will be 
detected, as a function of the central black hole spin assumed, and the assumed rate 
of EMR inspirals in M — 3 x lO 6 Af black holes. We show results for both detector 
configurations. The cases indicated by a * have z m i n > z max , but this is an artefact 
because in those cases less than 2 events are observed in total. 



in which phase the spin of the black hole plays a minor role. With pessimistic detector 
assumptions, the amount of observable inspiral is decreased, so the minimal redshift is 
larger, but still largely independent of assumptions on the black hole spin. 

It has been suggested that EMRI observations could be used to make high 
precision measurements of the Hubble constant in a statistical way (20] . However, these 
measurements would require LISA to detect ~ 10 events at z < 0.23 or ~ 20 events 




Figure 3. Distribution of LISA EMRI events as a function of the intrinsic central 
black hole mass (upper plot) and as a function of rcdshift (lower plot), for an inspiraling 
black hole with \i = 10Af Q . We show distributions for both detector configurations, 
computed using both the raw data and the fitting functions given earlier, as labelled. 



at z < 0.5. For the EMRI population we are considering here, assuming all MBHs 
have spin a = 0.9, we find that a rate 7^w > 14Gyr _1 will produce > 10 events at 
z < 0.23, and T^w > lOGyr" 1 will produce > 20 events at z < 0.5. If we assume 
all MBHs are non-spinning, a = 0, these increase to 15Gyr _1 and llGyr -1 respectively. 
The corresponding rates for the pessimistic detector configuration are 43Gyr~ 1 /41Gyr~ 1 
for a = 0.9 and 48Gyr _1 /56Gyr _1 for a = 0. Taking /?bh = rather than /?bh = —0.15 
these numbers become 21Gyr~ 1 /14Gyr~ 1 (a = 0.9) and 24Gyr~ 1 /16Gyr~ 1 (a = 0) or 
61Gyr~ 1 /54Gyr _1 (a = 0.9) and 73Gyr _1 /83Gyr _1 (a = 0) for the pessimistic detector 
configuration. These rates are all smaller than the current best estimates in [19], which 
indicates that this measurement should be feasible, although we start to face problems 
if the intrinsic rate is a factor of 5 lower and we have the pessimistic detector. 
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It has also been suggested that WD EMRIs could be used to make Hubble 
constant measurements, if an electromagnetic counterpart is observed from the tidal 
disruption of the WD [18J. This requires us to observe a WD inspiraling into a black 
hole with M < 10 5 M Q . Our prescription predicts only 0.1 WD events in the range 
1O 4 M < M < 10 5 M Q , which means it is unlikely that we would see an event. However, 
given the huge astrophysical uncertainties, we should not rule these out completely. 

4. Using LISA as an astrophysical probe 

For the simple model described in the previous section, we have seen that LISA should 
detect many EMRI events out to redshift z ~ 1.5. Each LISA event should provide 
us with the system parameters to high precision [3], and hence we will determine the 
masses and spins of as many as 1000 MBHs in the nearby Universe. Even if the intrinsic 
inspiral rate per black hole is much less or the LISA detector performs sub-optimally, we 
should characterize a few tens of low redshift MBHs. This information can be used for 
astrophysics. Due to the uncertainties in the rate of EMRIs per galaxy, the total number 
of events is not a good probe of the MBH population, but the distribution of events as a 
function of mass and/or spin is a useful probe, modulo the a + /3 7 degeneracy discussed 
earlier. It is for this application that the research described here was started and where 
the interesting astrophysics lies. In order to quantify what LISA can do, we are using 
the observable lifetime functions described here, in conjunction with results from black 
hole merger trees [21] to investigate the dependence of the LISA event distribution on 
the MBH population model. Crudely, if we have as many as 1000 events, then we can 
imagine dividing the mass range into ten bins, such that we would expect 100 events in 
each bin with our baseline model. The Poisson error in each bin would be ~ 10, and 
so if the Universe differed from our model in a way that changed the number of events 
in a bin by more than 10 we might hope to detect that in our observations. We are 
currently comparing different galaxy evolution models in a Bayesian way, by evaluating 
the likelihood that a given observed population of EMRI events came from a Universe 
with a particular history. This will allow an estimation of the number of EMRI events 
that will need to be observed in order to make a strong statistical statement about 
which model best describes the real Universe. This work is in progress, but we hope 
that the observable lifetime fits that we have presented in this paper will allow others to 
follow a similar research programme without the overhead of doing the LISA response 
calculations. 

Our present results have some limitations and we plan to quantify these in the 
future. One important question is how many of these events are "useful" in the sense that 
they give good parameter estimation. Naively we would hope that any event detected 
would have reasonable parameter estimation for the intrinsic parameters, but some 
nearby systems might be observed several years before plunge and so might not show 
enough evolution over our observation to allow accurate parameter determination. The 
fraction of the events that are useful should be quantified by a Fisher Matrix analysis. 
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Another issue is that our waveform model is restricted since we consider only circular- 
equatorial orbits. It is necessary to check how representative these are by considering 
observable lifetimes for eccentric and inclined orbits. As perturbative waveforms for 
generic orbits are computationally very expensive, a first cut at this problem would 
involve computing these lifetimes using approximate waveforms [3J E]. In addition, in 
the present work we are computing and using the observable lifetime of a sky-averaged 
source, but it should be verified explicitly that this does not differ too greatly from the 
sky-averaged observable lifetime of a source, which is what is really required. Finally, 
these SNRs are based on the low-frequency approximation to the LISA response and 
hence will not be completely accurate for systems with particularly low mass or high 
spin central black holes. This error needs to be quantified by computing more accurate 
SNRs for those systems using the full LISA response, although we do expect our current 
results to be conservative in this regard. 
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